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Abstract 

For the coherently driven A-type three-level systems the general ready- 
to-calculate expression for the susceptibility tensor at the frequency of the 
weak probe field is obtained for the arbitrary polarization of the strong 
coupling laser field and arbitrary values of the angular momenta of reso- 
nant atomic levels. The dependence of the relative difference in the group 
velocities of the two polarization components of the probe field on the 
polarization of the coupling field is studied. 

1 Introduction 

The phenomenon of the electromagnetically induced transparency (EIT), when 
the optical properties of the media for the weak probe field are coherently con- 
trolled by the strong coupling field, has been extensively studied in the recent 
years and provided a number of attractive applications [1 j . The most prominent 
feature of EIT is a spectacular reduction of the group velocity of light pulses 
in the resonant media [H [3] ■ The strong coupling field with the definite polar- 
ization produces also the optical anisotropy for the probe field, which reveals 
itself in the electromagnetically induced birefringence and polarization rotation 
of the probe field [4] [5]. Among the most promising applications of EIT is 
the implementation of quantum memory [HI El IB] - Such memory based on the 
controlled adiabatic deceleration and acceleration of single-photon pulses in the 
resonant media was suggested in [TO] , and soon the possibility of storage of 
light pulses was realized in rubidium vapor in the experiment The recent 
experiments [TS] [T3j [14] on EIT-based quantum memory demonstrate the con- 
tinuously improving efficiency and fidelity. There are different ways to encode 
the single-photon qubit, for example, in two spectral components of the photon 
pulse, as it was recently proposed in [15] , but the most natural way for qubit en- 
coding is provided by the photon two polarization degrees of freedom, as it was 
implemented in the experiments [T2l [13] . To store the photon polarization qubit 
its both polarization components must be effectively stopped in the medium, 
however the group velocities of these two components may differ essentially due 
to the optical anisotropy induced by the polarization of the coupling field. The 
experiments on EIT are performed on the the three-level A-type systems with 
degenerate levels, which are in many cases the hyperfinc structure components 
of alkali atoms degenerate in the projections of the atomic total angular mo- 
mentum on the quantization axis. The theoretical treatment of EIT on such 
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Figure 1: The level diagram. 



levels involves the solution of the equations for the atomic density matrix and 
the number of elements of this matrix increases drastically with the increase of 
the level angular momentum values. The objective of the present article is to 
obtain the general expression for the linear susceptibility tensor at the frequency 
of the probe field, which describes the polarization properties of EIT, and to 
study the dependence of the relative difference in the group velocities of the two 
polarization components of the probe field on the polarization of the coupling 
field. The essential point about this expression is that the susceptibility tensor 
may be easily calculated by means of standard matrix operations for the arbi- 
trary polarization of the driving field and for the arbitrary values of the level 
angular momenta. 



2 Basic equations and relations 

We consider the weak probe field propagating along the sample axis Z with the 
carrier frequency u>, which is in resonance with the frequency wo of an optically 
allowed transition J a — > J c between the ground state J a and the excited state 
J c , and the strong coherent coupling field propagating in the same direction 
with the carrier frequency lu c , which is in resonance with the frequency u> c o of 
an optically allowed transition Jb — >■ J c between the long-lived state Jb and the 
same excited state J c (Fig.l). Here J a , Jb and J c are the values of the angular 
momenta of the levels. The electric field strength of the coupling field (inside 
the sample as well as outside) and that of the probe field incident on the sample 
border z = may be put down as follows: 

E c = e c l c e-^( t - z / c )+c.c, (1) 

Eo=e e-^*-*/ c) +c.c., (2) 
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where e c and l c are the constant amplitude and the unit polarization vector 
of the coupling field, while eo is the constant vector amplitude of the incident 
probe field. At some sample point z, where the vector amplitude of the probe 
field is e, the evolution of the atomic slowly-varying density matrix p in the 
rotating-wave approximation is described by the equation: 



dp 



V,p 



(3) 



rel 

V = 2(AP c + SP b ) + n c (g c + gl) + G + G\ (4) 
g c = eX, G=(2\d\/h)ge*. (5) 

Here A = uj — ujq and A c = lo c — uj c o are the frequency detunings from resonance 
of the probe and of the coupling fields, while S = A — A c , P a is the projector 
on the subspace of the atomic level J a (a = a, b, c), £l c = 2\d c \e c /h is the 
reduced Rabi frequency for the coupling field, d — d(J a J c ) and d c = d(JbJ c ) 
being the reduced matrix elements of the electric dipole moment operator for 
the transitions J a — > J c and — > J Cl while g and g c are the dimensionless 
electric dipole moment operators for the transitions J a — > J c and Jf, — > J c . 
The matrix elements of the circular components g q and g cq (q — 0, ±1) of these 
vector operators are expressed through Wigner 3J-symbols 

(x amc = (-iy«--«(_i a I i), (6) 

\bc _ f _-i\Jb~m b I Jb 1 Jr. 



Finally, the term (dp/dt) re i in the equation ([3]) describes the irreversible relax- 
ation. Initially the atoms are at the ground state a the atomic density matrix 
being 

P(0) = Pa- 

Generally the atoms are at the equilibrium ground state with equally populated 
Zeeman sublevels, however they may be prepared at some special state [13] . 
In the linear approximation for the probe field we obtain from the equations 
([3]) for the elements of the steady-state (dp/dt = 0) atomic density matrix 
the following relations: 

( 7 -iA)p ca = l -ft c gtp ba + &p a ), (8) 
{T-i5)p ba = l -n c g c p ca , (9) 

where 

p a P = p a pp j3i a,P = a,b,c, 
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while the irreversible relaxation is simply characterized by the two real relax- 
ation rates - 7 for the optically allowed transition J a — > J c and T for the 
optically forbidden transition J a — >• J5: 

ia\ co / 7 ~ \ ba 

i) (I) ,- r ^ 

/ rel \ / rei 

From (J8j)-([9| it follows immediately: 

r - sj^e-^A., do) 

where 

^ = ^ + 4( 7 -4r-^) ^ c - (U) 

The medium polarization component at point z 

P = pe-' lu}t + c.c. 

with the frequency uj of the probe held is expressed through the density matrix 
p ca defined by (fT0 l) -(fTT ]) : 

p = n \d\Tr A {gp ca }, (12) 

where Uq is the concentration of resonant atoms and the trace is carried out over 
atomic states. With the two orthonormal vectors 1^ in the polarization plane 
XY (Ijll = S jk , j, k = 1, 2) the equation (O with an account of (T0 ]) -(fTT |) and 
([5]) may be expressed through the susceptibility tensor Xjk- 



Pj 



= £ o^2xjke k , (13) 



fc=i 

where 

X]k = \-7{kh) TrA {^ _1 $*a>&} ' (i4) 

n o\ d \ 2 * „,* , 1K x 
Xo = — r— , & = glj- 15 

With the introduction of the orthonormal set of eigenvectors \c n > with non- 
negative eigenvalues (? n of the hermitian operator g\g c '- 

g\g c \cn >= cl\c n >, n = 1, 2 J c + 1, 

the susceptibility tensor (fT4|) may be transformed to the expression: 

*xo < c n \g' k Pagj\c n > 



Xjk 



1 - i(A/ 7 ) £j 
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c 2 n 2 

A " = 1 + 4( 7 -zA)(r^r (17) 

Now let Xk be the two eigenvalues and 

2 

vfc=5^» 3 -fclj, fc = l,2, (18) 

3=1 

be the corresponding two eigenvectors of the 2x2 susceptibility matrix x = 
{Xjk}- 

XVfc = XfcVfc, k = 1,2. 

Then, the vector amplitude eo of the incident probe field after the passage of 
the distance z in the medium is transformed to 

e = 5e , S = vfv~\ (19) 

where v = {vj k } is the 2 x 2-matrix, defined by the equation (TT5|) . while the 
diagonal matrix T is determined by the eigenvalues Xk of the susceptibility 
matrix: 

T jk = 5 jk e^' c ^ , nj = y/T+xl, (20) 
as it may be obtained in a usual way from the Maxwell equation 

<9 2 E 1 <9 2 E _ 1 d 2 P 

dz 2 c 2 dt 2 eqc 2 dt 2 

for the electric field strength E of the probe field. The intensity I and the 2x2 
polarization matrix a — {&jk} of the probe field after the passage of the distance 
z in the medium are related to the intensity Iq and the polarization matrix ctq 
of the incident probe field by the equations: 

f =tr(Sa &),a= (21) 



where transformation matrix S is defined by (|19 p -(|20 |) . 

The electromagnetically induced transparency is associated with the signifi- 
cant reduction of the group velocity of the probe field due to the steep dispersion 
at the transparency window. Since the dispersion is different for the two polar- 
ization components of the probe field, these two components will be differently 
slowed down in the medium. Let us calculate the relative difference of the 
group velocities for these two polarization components. For the observation of 
the electromagnetically induced transparency the coupling field must be strong 
enough: 

n 2 c » r 7 . 

Then, close to the two-photon resonance 5 < T the general formula (ITBl) - (|17p 
for the susceptibility tensor is simplified to the following expression: 
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Figure 2: The imaginary part of susceptibility Im(xxx/xo) (clashed line) and 
Im(xyy/xo) (solid line) versus dimensionless frequency detuning A/7 on the 
transitions with J a — 0, J& = J c = 1 at f2 c /7 = 1, r/7 = 0.01 and l c = l x . 
A c = 0. 



where 

< c n \g\p a g 3 \c n > 
a jk = -2 . ( 23 ) 

n n 

is a hermitian 2x2 matrix. The summation in (|23[) is carried out over all 
eigenvectors \c n > of the operator g\g c with non-zero eigenvalues c\ > 0. The 
two eigenvectors (jT5)) of the susceptibility matrix {Xjk} coincide with the two 
orthonormal eigenvectors of the hermitian matrix {ajk}, while its eigenvalues 



Xk = -^2~ {5 + iT)a k , 



are expressed through the real eigenvalues a k of matrix {a.,*,}. The group ve- 
locities V? of the two probe field components polarized along eigenvectors 
are determined by the refraction indices n' k = Re (-^1 + Xk) '■ 



n' k + uj(dn' k l 'duj) 
Since \xk\ <C 1, in the case of steep dispersion 

Xoiuj > tt 2 c , 
we obtain for the group velocities: 



V*-— r. 9r - 2xo ^ ak (04) 



The relative difference of these two group velocities 

vr - vr 
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Figure 3: The real part of susceptibility Re(xxx/xo) (dashed line) and 
Re(x vy /xo) (solid line) versus dimensionless frequency detuning A/7 on the 
transitions with J a — 0, J b = J c = 1 at £l c /l = 1, r/7 = 0.01 and l c = \ x , 
A c = 0. 

is determined by the difference of the two eigenvalues of matrix {a^fc}, ai be- 
ing the largest of them and the group velocity Vf of the probe field component 
polarized along the corresponding eigenvector Vi being the smallest of two. 

3 Discussion 

The (2J b + 1) x (2J C + 1) matrices g c and (2J a + 1) x (2J C + 1) matrices gj 
(j = 1,2) of the electric dipole moment operators may be easily calculated for 
any reasonable values of the angular momenta J a , Jb and J c of the resonant 
atomic levels using the known formulae for Wigner 3J-symbols, and further cal- 
culations of the susceptibility tensor x (|14[) and of the transformation matrix S 
P9[) . involving standard matrix operations, are rather simple. In the following 
numerical calculations we shall assume the coupling field to be rather strong: 
fl c /j — 1 and exactly resonant A c = 0, and the relaxation rate of the forbidden 
transition to be rather small: r/7 = 0.01, as it is the case for the electromag- 
nctically induced transparency, while the atoms are initially at the equilibrium 
state: 

. Pa 

Pa ~2J a + l- 

We shall also consider the linearly polarized coupling field, then without loss of 
generality its polarization vector may be directed along the Cartesian axis X: l c 
= la:, and we shall choose the Cartesian basis li = l x , l2=ly in the polarization 
plane XY for the calculation of the elements of the susceptibility tensor. With 
this choice of the Cartesian basis the polarization matrix a (f2"Tj) of the probe 
field may be expressed through the Stokes parameters (n — 1, 2, 3): 

„ = l( 1 + 6 

2 Ui+to 1-6 /' 



7 




medium depth z/L 



Figure 4: The relative intensity I/Iq (dashed line) and the degree of polariza- 
tion P (solid line) of the initially unpolarized probe held versus dimensionless 
medium depth z/L (L — l/kxo) on the transitions with J a = 0, = J c = 1 at 
ft c / 7 = 1, r/ 7 = 0.01 and l c = l x , A = A c = 0. 



p = ^l + Q + ti 

being the total degree of polarization, 
being the degree of linear polarization, and 

Pc=U 

being the degree of circular polarization. 

Let us consider the transitions with the low values of angular momenta 
J a = 0, Jb = J c = 1, which may be realized in thallium vapor for example. 
With the linear polarization l c = 1^ of the coupling field the susceptibility ten- 
sor becomes diagonal in the Cartesian basis. The dependencies of the imaginary 
Im(Xxx/xo), Im(Xyy/xo) and real Re(x X x/xo), Pe(x yy /xo) parts of the suscep- 
tibility eigenvalues versus the dimensionless frequency detuning A/7, obtained 
from the calculations according to the formula (fl"6|) . are presented in the Fig- 
ures 2 and 3. As it follows from these dependencies the strong linearly polarized 
(along the axis X) coupling field opens the "transparency window" (solid line) 
for the polarization component of the probe field linearly polarized in the per- 
pendicular direction (along the axis Y), while for the component of the probe 
field collinearly polarized with the coupling field (along the axis X) this "win- 
dow" remains shut (dashed line), which gives rise to rather strong dichroism 
for the resonant (A = 0) probe field. This happens because the X-polarized 
component of the probe field interacts with a single substate 



\c x >= — (\m c = -1 > -\m c = 1 >) 
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Figure 5: The imaginary part of susceptibility Im(Xxx/xo) (dashed line) and 
Ifn(xyy/Xo) (solid line) versus dimensionless frequency detuning A/7 on the 
transitions with J a = 1, Jf, = J c = 2 at £l c /j = 1, r/7 = 0.01 and l c = L,., 
A c = 0. 



of the excited level c, which is not affected by the coupling field: 

g c \c x >= 0. 
For the initially unpolarized probe field: 

the dependencies of its relative intensity I/Iq (dashed line) and the degree of 
polarization P (solid line) versus the dimensionless medium depth z/L, where 
L = 1/kxoi are presented in the Figure 4. As it may be seen in this figure, 
the probe field becomes almost fully polarized with the degree of polarization 
P = 0.96 (linearly along the axis Y), while its intensity is reduced by the factor 
of 0.19. 

Let us now consider the transitions with larger values of angular momenta 
J a = f, Jh = J c = 2, which were employed in the experiment [4] performed 
in 87 Rb vapor. For the linear polarization l c = \ x of the coupling field the 
dependencies of the imaginary Im(x xx /xo), Im{Xyy/Xo) and real Re(x xx /xo), 
Re{Xyy/Xo) parts of the susceptibility eigenvalues versus the dimensionless fre- 
quency detuning A/7 are presented in the Figures 5 and 6. In this case the 
"transparency windows" are open for both polarization components of the probe 
field, however, that for the ^-component (solid line), perpendicular to the cou- 
pling field, is larger than for the X-component (dashed line), collincar with the 
coupling field. In the case of large values of the angular momenta the dichroism 
is less than in the case of small ones, as it may be seen in the Figure 7. When 
the degree of polarization of the initially unpolarized probe field attains the 
value of 0.5 its intensity is reduced by the factor of 0.17. 
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Figure 6: The real part of susceptibility Re(xxx/xo) (dashed line) and 
R e iXyy/Xo) (solid line) versus dimensionless frequency detuning A/7 on the 
transitions with J a = 1, Jf, = J c = 2 at £l c /j = 1, r/7 = 0.01 and l c = L,., 
A c = 0. 

The relative difference e in the group velocities of the two polarization com- 
ponents of the probe field (l2"5j) depends only on the values of total angular 
momenta of resonant levels, on the initial atomic state and on the polarization 
of the coupling field. Let us consider the coupling field with arbitrary elliptic 
polarization in the plane XY: 



where parameter a defines the ratio of the lengthes a y and a x of the axes of 
polarization ellipse along the axes Y and X according to the relation 



while the sign of tan(a — 7r/4) determines the direction of rotation of the pulse 
electric field vector in the plane XY. For the atoms at initially equilibrium state 
and with the fixed values of the level angular momenta the relative difference in 
the group velocities e depends only on the ellipticity parameter a of the coupling- 
field. An example of such dependence for the levels with the angular momenta 
Ja = Jc = 1) Jb = 2, employed in the experiment [T3], is presented in the Figure 
8. The maximum difference in the group velocities e = 0.625 is obtained with the 
circularly polarized coupling field (a = 0,7r/2), while the minimum difference 
e = 0.125 is obtained with the linearly polarized coupling field (a — tt/4). 
In the case of circularly polarized coupling field the polarization component 
with the larger group velocity is circularly polarized in the same direction as 
the coupling field, while the polarization component with the smaller group 
velocity is circularly polarized in the direction opposite to that of the coupling 
field. In the case of linearly polarized coupling field the polarization component 
with the smaller group velocity is linearly polarized in same direction as the 
coupling field, while the polarization component with the larger group velocity 
is linearly polarized in the direction perpendicular to that of the coupling field. 
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Figure 7: The relative intensity I/Iq (dashed line) and the degree of polariza- 
tion P (solid line) of the initially unpolarized probe field versus dimensionless 
medium depth z/L (L = \/kxo) on the transitions with J a = 1, = J c = 2 at 
r> c / 7 = 1, r/ 7 - 0.01 and l c = \ x , A = A c = 0. 

In the experiment [13] the coupling field was linearly polarized in the direction 
of propagation of the probe field: l cq = <5 g .o (7r-polarized) , while it propagated in 
the perpendicular direction, and the atoms were prepared at the pure Zeeman 
state | J a — 1, rn a — > with the zero projection on the quantization axis. Under 
such conditions we obtain for the tensor (|23|) : ajk = 1.667(5^, so that the group 
velocity of the probe pulse in this case does not depend on its polarization. For 
the unprepared atoms, which are initially in the equilibrium state with equally 
populated Zeeman sublevels, we obtain ajk = 0.972djk, so that in this case the 
group velocity of the probe pulse also does not depend on its polarization. 

4 Conclusions 

In the present article the A-type three-level systems with degenerate levels coher- 
ently driven by the strong laser field are considered. The general expression for 
the susceptibility tensor at the frequency of the weak probe field is obtained for 
the arbitrary polarization of the strong coupling laser field and arbitrary values 
of the momenta J a , J^, J c of resonant atomic levels. The numerical calculations, 
based on this expression, may be easily performed by means of standard matrix 
operations for any reasonable values of angular momenta. Sample calculations 
for the transitions with J a — 0, J& = J c = 1 and J a = 1, Jf, = J c = 2, which 
are employed in the experiments with atomic vapors, were carried out. It was 
shown, that the depth of the "transparency window" for the resonant probe 
field depends on its polarization providing rather strong dichroism. Such de- 
pendence is greater for the low values of angular momenta J a — 0, Jf, = J c = 1 
and it is muffled for larger values J a = 1, J& = J c = 2. The general expression 
for the relative difference in the group velocities of the two polarization com- 
ponents of the probe field is obtained. This difference depends on the values of 
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Figure 8: The dependence of relative difference e = (V^ — Vf r )/V~2 r in the 
group velocities of the two polarization components of the probe field on the 
ellipticity parameter a of the coupling field on the transitions with J a = J c = 1, 
Jb = 2. 

total angular momenta of resonant levels, on the initial atomic state and on the 
polarization of the coupling field. Such dependence on the ellipticity parameter 
of the coupling field for the transitions with J a = J c — 1, Jb — 2, studied nu- 
merically, revealed rather a wide range of variation of this relative difference - 
from 0.625 with the circularly polarized coupling field to 0.125 with the linearly 
polarized one. 
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